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The standard description of Fermi acceleration, developing in a class of time-dependent billiards, is given in 
terms of a diffusion process taking place in momentum space. Within this framework the evolution of the 
probability density function (PDF) of the magnitude of particle velocities as a function of the number of 
collisions n is determined by the Fokker-Planck equation (FPE) . In the literature the FPE is constructed by 
identifying the transport coefficients with the ensemble averages of the change of the magnitude of particle 
velocity and its square in the course of one collision. Although this treatment leads to the correct solution after 
a sufficiently large number of collisions has been reached, the transient part of the evolution of the PDF is not 
described. Moreover, in the case of the Fermi-Ulam model (FUM), if a stadanrd simplification is employed, 
the solution of the FPE is even inconsistent with the values of the transport coefficients used for its derivation. 
The goal of our work is to provide a self-consistent methodology for the treatment of Fermi acceleration in 
time-dependent billiards. The proposed approach obviates any assumptions for the continuity of the random 
process and the existence of the limits formally defining the transport coefficients of the FPE. Specifically, we 
suggest, instead of the calculation of ensemble averages, the derivation of the one-step transition probability 
function and the use of the Chapman-Kolmogorov forward equation. This approach is generic and can be 
applied to any time-dependent billiard for the treatment of Fermi-acceleration. As a first step, we apply this 
methodology to the FUM, being the archetype of time-dependent billiards to exhibit Fermi acceleration. 
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Fermi acceleration, that is the increase of the 
mean energy of an ensemble of particles due to 
random collisions off moving scatterers, is clearly 
one of the most interesting physical mechanisms 
linked to time-dependent billiards. Despite this 
fact, the standard approach in the literature 
for its analytical treatment can at best describe 
Fermi acceleration in the asymptotic time limit. 
Herein, we propose a methodology, which de- 
scribes the evolution of Fermi acceleration at all 
times and, even more, obviates any unclear or ad 
hoc assumptions, which can lead to inconstencies 
or unreliable results. We exemplify the proposed 
approach in the prototype of billiards exhibiting 
Fermi acceleration; The Fermi-Ulam model. 



I. INTRODUCTION 

More than 60 years ago, Fermii proposed an intuitive 
mechanism for the explanation of the origin of the highly 
energetic cosmic ray particles and ever since it has been 
a subject of intense study. The mechanism consists in 
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the increase of the mean energy of particles as a result 
of random collisions with moving scatterers. Soon after 
his seminal paper, his co-worker Ulam introduced a sim- 
ple mechanical model for testing Fermi's idca^, known as 
the Fermi-Ulam model (FUM), linking for the first time 
Fermi acceleration with the study of time-dependent bil- 
liards. 

Since the introduction of the FUM, the standard de- 
scription of Fermi acceleration developing in a class of 
time-dependent billiards is given in terms of a diffusion 
process taking place in momentum space^Ti^. Within this 
framework the evolution of the probability density func- 
tion (PDF) of the magnitude of particle velocities as a 
function of the number of collisions n is determined by 
the Fokker-Planck equation (FPE) . In the literature the 
FPE is constructed by identifying the transport coeffi- 
cients with the ensemble averages of the change of the 
magnitude of particle velocity and its square in the course 
of one collision^Ti^. Although this treatment leads to the 
correct solution after a sufficiently large number of colli- 
sions has been reached, the transient part of the evolution 
of the PDF is not described. Moreover, in the case of the 
FUM, if a standard simplification is employed — known 
as the static wall approximation (SWA) or the simplified 
Fermi Ulam Model (SFUM)— the solution of the FPE is 
even inconsistent with the values of the transport coeffi- 
cients used for its derivation. 

The aim of the work presented, is to provide a self- 
consistent methodology for the derivation of the PDF 
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of particle velocities for all times. The proposed ap- 
proach obviates any assumptions for the continuity of the 
random process and the existence of the limits formally 
defining the transport coefficients of the FPE. Specifi- 
cally, we suggest, instead of the calculation of ensem- 
ble averages, the derivation of the one-step transition 
probability function (TPF) and the use of the Chapman- 
Kolmogorov (forward) equation (CKE). This approach 
is generic and can be applied to any time-dependent bil- 
liard for the treatment of Fermi-acceleration. As a first 
step, we apply this methodology to the FUM, being the 
archetype of time-dependent billiards to exhibit Fermi 
acceleration. In this context, we show that the FPE re- 
ported in the literature^ describing the evolution of the 
PDF of the magnitude of particle velocities is not valid, 
and that the observed agreement for n ^ 1 between the 
analytical and numerical results, in this case, should be 
regarded as accidental, i.e. due to the validity of the 
central limit theorem (CLT). 



II. STATISTICAL DESCRIPTION OF FERMI 
ACCELERATION 



A. The Fokker-Planck approximation 

The standard approach in the literature for the de- 
termination of the asymptotic behaviour of the PDF of 
particle velocities, is the approximation of the jump pro- 
cess with a diffusion process, described by the FPE^"— . 
This approximation has also been used for the analyti- 
cal treatment of Fermi acceleration developing in higher- 
dimensional billiards, like the simplified periodic Lorentz 
gag^, i.e. the oscillating hard circular scatterers oscillate 
only in the velocity space. An equation of the form of 
the FPE can be derived from Eq. ([T} as follows'': 

If we introduce Av = w — z, then the integrand in 
Eq. ID) can be rewritten as. 



(2) 



p {v,n\v' ,n') = J d{Av)p {v — Av,n — l\v' ,n') 
X W {v — Av + Av, V — Av) . 



Expanding the distribution function p (v, Av,v' ,n') and 
the transition probability function (TPF) W(v; Av) in a 
Taylor series yields. 



Fermi acceleration developing in a time-dependent bil- 
liard can be described in terms of a stochastic process 
taking place in the velocity space. Let W{v, z) denote 
the probability of a particle being at the velocity z to 
perform a jump to velocity v in the course of a single 
collision and p{v, n\v' , n') the probability of a particle to 
possess velocity v after n collisions given that at the n'th 
collision it had velocity v' . This jump process can be 
described by the following equation: 

p {v,n\v' ,n') — / dzp {z,n — l\v' ,n')W {v, z) (1) 



Eq. ID) is exact, on the condition that the process is 
Markovian. From a physical point of view, this means 
that the probability of a particle to experience a velocity 
jump equal to Av upon the nth collision depends only on 
the velocity it had at the previous step, i.e at the nth-1 
collision. 



r oo / 1 \m 

p{v,n\v\n') = / diAv) ^ ^-^{Av) 



X Q-^P n - l\v', n) W {v + Av, v) 



(3) 



Integrating now Eq. |2|) over Av we obtain, 



p{v,n\v',n') = ^ 



(-1)™ d" 



m— 



i! dv' 



-Mm{v)p{v,n- l\v' ,7l') , 
(4) 

where Mm{v) stands for the mth moment of the TPF, 
i.e. 



Mm{v) = / iAv)"'W{v + Av, z)d{Av). 



Therefore, 



p {v, n\v' , n) — p{v,n ~ l|i>', n) 



(-1)™ d" 
ml dv^ 



:M„,{v)p{v,n\v',n') 



(5) 



By truncating the above series to the second order, and for n ^ 1 one obtains an equation resembling the FPE. 
further by approximating the discrete derivative 

d d 
—p {v, n\v', n')=^- — [Bp [v, n\v' , n')] 
an ov 

Akp{v,n\v',n')^[p{v,n + k\v',n')- p{v,n\v',n')]lk, I d'' , ^ , ^' 



(fc = 1) with the continuous derivative dp{v, n\v' , n') / dn, where the coefficient B and D is the ensemble average of 
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the change of particle velocities and its square, respec- 
tively, in one mapping period. 

The approximations applied above for the construction 
of the FPE are valid on the condition that only very 
small jumps are probable and further that the solution 
p (u, n\v' ^ n') varies slowly with v so that one can perform 
the expansion in a Taylor series. More formally^, we 
demand that there exists a (5 > 0, 

T4^(z + Az, z) « 0, for |Az| > (5 (7a) 

p{v + Av, n\v', n') « p{v, n\v\ n'), for | Au| < S. (7b) 

In the literature^, the derivation of an FPE from Eq. ([1]) 
for the statistical description of Fermi acceleration is car- 
ried out on an ad hoc basis. As a consequence, as shown 
in the following, it has produced contradictory results. 
Moreover, by construction, the description of Fermi ac- 
celeration with a continuous stochastic process, can at 
best describe the statistics only for n ^ 1. Hence, a full 
description of FA in a time-dependent billiard can only be 
given in the context of a jump process and consequently 
by Eq. ©. 



Moreover, in this case an approximate solution can be 
obtained directly in the velocity space, using the sad- 
dlepoint approximation technique^. Specifically, from 
Eq. ([5]), one can derive the moment generating function 



e*''^p{v, n\z, 0)dv 



e'''W{v)dv 



(11) 



of the velocity PDF. To find the saddlepoint t{v,n), 
we solve the equation K!{t,n) = v, where K{t.n) — 
\og{(j){t,n)). Then, the PDF is approximately. 



p{v,n\z,0) 



y 2TrK"{t{v,n)) (12) 
X exp [k(£(u, n)) — i{v, n)u] . 



In the following sections we implement the proposed 
methodology in the prototype of time-dependent billiards 
exhibiting Fermi acceleration; The Fermi-Ulam model 
(FUM). 



B. A complete description: The Chapman- Kolmogorov 
equation 

The study of the transient statistics can only be accom- 
plished by means of the Chapman-Kolmogorov equation, 
i.e. Eq. ([Ij. Assuming that initially particle velocities are 
distributed according to p{v,0) — 6{v ~ z), Eq. ([T]) can 
be rewritten in respect with the one-step TPF. W{v, v') 
as, 

p(?;„,n|z,0) = j ■■■ j W{vn,Vn-i)---W{vi,z)dw, (8) 

n-l 

where dv = Y\ dvi. The derivation of the one-step 

TPF can be achieved by determining the PDF p{q) of 
the variables q = {xi} appearing in the dynamical equa- 
tion defining the velocity of a particle after a collision 
with the moving boundary of the time-dependent billiard, 
Vn = f{vn-i,q)- Then, the TPF is 

W{Vn,Vn^l) = J p{q)6[vn - fiq,Vn-i)]dq. (9) 

If the resulting TPF is a function of the difference of ve- 
locities at successive steps M^(w„,w„_i) — W{vn — Vn-i), 
Eq. ^ can be easily solved in the Fourier space. Specif- 
ically, if this condition is met, then by taking the Fourier 
transform of Eq. ([S]) we find, 

J- [p{v, n\z, 0)] = (27r)^ e'^'^^ [W{v)]r , (10) 

OO 

where = l/(-\/27r) J exp[~-ikv]dv. 

— CX3 



III. FERMI ACCELERATION IN THE STOCHASTIC 
SIMPLIFIED FUM 

The Fermi-Ulam model, originally proposed for test- 
ing the feasibility of gaining energy through scatter- 
ing off moving targets, i.e. Fermi acceleration, consists 
of one harmonically oscillating and one fixed infinitely 
heavy hard wall and an ensemble of non-interacting par- 
ticles bouncing between them. Ever since, many dif- 
ferent versions of the original model have been sug- 
gested and investigated, such as variants of the FUM 
with dissipationi^^ii^ different deterministic or random 
drivings of the moving wal l^^i^^ the quantum-mechanical 
versionii""S^ and the, so called, bouncer model^^, where 
a particle performs elastio^i or inelastio^"— collisions 
with an oscillating infinitely heavy platform under the 
influence of a gravitational field. Recently, a hybrid ver- 
sion of the FUM and the bouncer model has also been 
investigated^^i^. 

The equations defining the dynamics of the FUM are 
of implicit form with respect to the collision time, which 
complicates numerical simulations and hinders an analyt- 
ical treatment. A simplification^ — known as the static 
wall approximation (SWA) -'^^i^^ consists in treating the 
oscillating wall as immobile, located at its equilibrium 
position, yet allowing the transfer of momentum upon 
impact with a particle as if the wall were harmonically 
oscillating. This simplification has become over the time 
the standard approximation for studying the FUM'^-. 
The SWA speeds-up numerical simulations and facili- 
tates the analytical treatment of the problem, while it 
has been generalized to higher-dimensional billiards with 
time-dependent boundaries, such as the time-dependent 
Lorentz Gaa^^. 
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Let us consider, without loss of generality, a FUM con- 
sisting of a fixed wall on the right and a moving wall 
on the left, oscillating with frequency w. If we further 
assume that the positive direction of particle velocities 
is towards the right, then the dynamics of the billiard 
within the framework of the SWA is defined by the fol- 
lowing set of dimensionless difference equations: 



^n— 1 



Vn-1 
V„ = \Vn-l + 2Un\ 
Un = ecos(t„ -I- 77„), 



(13a) 

(13b) 
(13c) 



where u„ is the velocity of the "oscillating" wall, w„ is the 
algebraic value of the particle velocity immediately after 
the nth collision with the "oscillating" wall measured in 
units of LOW (w denoting the spacing between the walls), 
tn the time when the nth collision occurs measured in 
units of 1/w, r]n a random variable uniformly distributed 
in the interval [0, 27r) updated immediately after each 
collision between a particle and the fixed wall and e the 
dimensionless ratio of the amplitude of oscillation to the 
spacing between the "oscillating" and the fixed wall. It 
is noted that in all numerical simulations e was fixed at 
1/10. 

The absolute value in Eq. (jl3bp is introduced in order 
to avoid the occurrence of positive particle velocities after 
a collision with the "oscillating" wall, which would lead 
to the escape of the particle from the area between the 
walls. It should be stressed that such a collision, within 
the framework of the exact model, corresponds to a parti- 
cle experiencing at least one second consecutive collision 
with the "oscillating" wall. Therefore, if < 2|u„| 

and M„ < 0, in order to prevent the particle from escap- 
ing the region between the walls the velocity is reversed 
artificially. The presence of the absolute value function in 
Eq. (jl3b[) . nevertheless, complicates the analytical treat- 
ment of the acceleration problem. For this reason, it has 
become a standard practice in the treatment of the FUM 
to remove the absolute value function, thereby neglect- 
ing the set of collision events upon which the particle 
direction is not reversed after its collision with the "os- 
cillating" wall. However, this further simplification gives 
rise to a fundamental inconsistency: the ensemble mean 
of the absolute velocity obtained analytically does not 
change through collisions with the oscillating'^ wall, de- 
spite the well-established numerical result that Fermi ac- 
celeration does take place in the phase-randomized FUM. 



A. The asymptotics of the PDF of particle velocities 

1. Application of the central limit theorem 

In this section we will discuss the asymptotic behaviour 
of the PDF of particle velocities in the SFUM. Evidently, 
after n collisions the velocity of a particle evolving in the 
SFUM is the sum of the velocity jumps it has experi- 



enced up to this point, i.e. Vn = ^'^n + ''^o- Fur- 

m— 1 

thermore, due to Fermi-acceleration developing in the 
SFUM, after n ^ 1 collisions, the vast majority of the 
particles has acquired velocities much greater than the 
maximum wall velocity, irrespective of the initial distri- 
bution. Therefore, most of the collisions, after a suf- 
ficiently large "time", take place in the high velocity 
regime. In this limit, the absolute value function can be 
neglected and we immediately obtain^i^ i^^) = and 
((Av)^) = 2e^. Therefore, the velocity jumps are com- 
pletely uncorrelated, i.e. do not depend on the velocity 
the particle had at the previous step. Thus, for n ^ 1 
and V ^ e the central limit theorem (CLT) dictates that 
the PDF of particle velocities tends to a Gaussian distri- 



bution, with a mean value equal to ^ (Aw^ 

i=l 



and variance 



{Avi 



Hence, the PDF of particle velocities for n 1 is 



p{v,n) 



1 



cxp 



(14) 



In Fig. [U Eq. (fT4|) is plotted along with the histogram of 
particle velocities obtained from the simulation of 1.2 x 
10^ trajectories for n — 10^ collisions. The ensemble was 
initially distributed according to the delta function 6(v — 
e). The analytical result obtained from the application 
of the CLT is in perfect agreement with the numerically 
computed PDF. 



2. FPE equation in the SFUM 

As mentioned in the previous section, assuming that 
for n ^ 1 the probability measure of the events occurring 
in the low-velocity regime is negligible, B = (Av) i2 0, 
D = {{Av) ) i2 2e^. In this limit Eq. ^ obtains the form 
of a standard diffusion equation, which for a delta initial 
distribution of velocities v = z together with refiecting 
conditions at v = has as a solution the sum of two 
spreading Gaussians 



p{v,n\z,0) 



1 



2V 



TTne^ 



exp 



-I- exp 



{v - zf 



{v + z) 
4ne^ 



(15) 



which for n > z^/(4e'^ln2) transforms to Eq. (HH). 



3. Remarks 

Although the solution derived by means of the FPE 
is in agreement with the one obtained from the applica- 
tion of the CLT, the methodology used for the derivation 
of Eq. ^ stands on very shaky ground, since the ter- 
mination of the series at the second term in Eq. ([SJ is 
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FIG. 1. Histogram — diagonal crosses — of particle velocities 
after n = 10^ collisions, obtained by the iteration of Eqs. (|13|l . 
on the basis of an ensemble of 1.2 x 10^ particles initially 
distributed as p{v,0) = S{v — e). The analytical result derived 
through the application of the CLT [Eq. (|14p ] is also plotted 
— solid (red) line. 



completely arbitrary'^. In general, a jump process can 
be approximated by a diffusion process, on the condition 
that a scaling assumption for the transition probability 
holds. Namely, in the limit of infinitely small time in- 
tervals, the jumps should become smaller and more fre- 
quent, such that the random process can be viewed as a 
continuous one^. An intuitive way to examine this is to 
consider the average square of the jump size ((Av)^) a 
particle makes having a velocity v prior to the collision. 

Given that the SWA treats the moving wall as fixed 
in the configuration space, all phases upon collision are 
possible, independently of the velocity v of particles be- 
fore a collision. As a result, the average jump size is not 
reduced as u — >■ 0. In contrast, within the exact model, as 
the velocity of the particle prior to a collision decreases, 
it becomes increasingly probable to collide with the wall 
at the turning points, where its velocity is close to zero. 
Moreover, if the velocity of the particle before a colli- 



sion is small, then successive collisions are likely to oc- 
cur, the exact dynamics result to higher exit velocities. 
Consequently, successive collisions render small particle 
velocities improbable, as opposed to the SFUM, where as 
shown in Sec. lIII All u = is the most probable velocity. 

Summarizing, the application of the CLT for the de- 
termination of the long-time statistics is much more 
straightforward and renders the solution of a differential 
equation redundant. More importantly, the assumption 
of continuity of the stochastic process describing Fermi 
acceleration, which is essential for the construction of an 
FPE, is not required. 



B. Short-time statistics in the SFUM 

From Eqs. the particle velocity after the nth col- 
lision given that it had velocity z is. 



-z-2u-2{z + 2u) e (2u -I- z) 



(16) 



where Q{x) is the Heaviside unit-step function. 

According to Eq. (|13bp the wall velocity u„ is deter- 
mined by the phase ^„ = tn + rjn of oscillation at the 
instant of the nth collision. Due to the fact that in the 
stochastic SFUM the phase is randomly shifted through 
the addition of a random number tjn — distributed uni- 
formly in the interval (0, 27r) — after each collision, the 
oscillation phase ^„ is completely uncorrelated between 
collisions, following a uniform distribution. Furthermore, 
given that in the context of the SFUM the wall remains 
fixed in the configuration space, the wall velocity upon 
collision does not depend on the velocity of the particle, 
therefore, m„ and f„_i are also uncorrelated. From the 
fundamental transformation law of probabilities the PDF 
of the wall velocity upon collision is. 



p{u) = 



1 



(17) 



For the single-step TPF W{v, z) we can write, 

W{v,z)^J p{u)S[v-v{u,z)]du. (18) 

Substituting Eqs. (HZ]) and (HH) into Eq. ^ we obtain 
after integrating over u, 



W{v,z) 



e{2e-v - z,2e- z) 
^4e2 - (u + z)2 
e{2e-v + z,2e + 



^4e2 - (u - zf 



(19) 



In Fig. [2] the analytical result of Eq. (|19|) is compared 
with the histogram of particle velocities after a single 
collision, obtained numerically using Eqs. (|13l) and an 
ensemble of 1.2 x 10^ particles, with initial velocity z — 
0.1. Clearly, the numerical and analytical results are in 
agreement. 
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FIG. 2. Histogram — diagonal crosses — of particle velocities 
after a single collision with the "moving" wall, obtained using 
Eqs. (|13p and an ensemble of 1.2 x 10^ particles with initial 
velocity z = 0.1. The analytical result [Eq. (|19p ] for the 
one-step transition probability is also plotted for the sake of 
comparison — solid line 



The analytical result of Eq. ([H]) reveals that the TPF 
depends only on the most immediate history of a parti- 
cle, that is on the velocity it had at the previous step. 
Consequently, the stochastic process is indeed Marko- 
vian. Even more, if the particle before a collision has 
velocity z > 2e, then the velocity jump Aw — v ~ z it un- 
dergoes is completely independent on its history. There- 



fore, changes in velocity in the high-velocity regime are 
completely uncorrelated. 

In more detail, Eq. ([T9l) . consists of two parts, one of 
which does indeed depend only on the jump size. How- 
ever, the other branch of the TPF, taking effect for v < 2e 
— relating to the set of rare events^^ — depends also on 
the velocity at the last step. Nevertheless, the action of 
both branches of W allows of a simple geometrical inter- 
pretation: At each step, the second branch of the TPF 
stretches the PDF p{v,n\z,0), resulting to a probability 
flux towards negative values of velocity. This unphysi- 
cal result caused by the stretching is negated by the first 
branch, which folds the part of the p density over the 
vertical line at w = 0. Therefore, the solution of Eq. (jS]) 
can be obtained by extending the domain of p(v,n\z,0) 
to the whole real line and applying the method of images. 
Thus, for any number of collisions, we have 



p{v, n\z,0) ~ p{v, Ti\z,0) + p{v, n\ — z, 0), 



(20) 



where p is the solution of the unrestricted problem. Sub- 
stituting Eq. (ITOl) into Eq. we obtain. 



p(fc, n\z, 0) = exp {—ikz)Jo (2e|/c|)' 



(21) 



Eq. ([2T|l cannot be inverted analytically. To obtain an 
analytical result into the velocity space, we use the sad- 
dlepoint approximation [Eq. (jl2p ]. The moment gener- 
ating function of p{v, n\z, 0) is, 

0(^,^^) = (/o(2^e))"e*^ (22) 

where Iq is the modified Bessel function of the first kind. 
Consequently, the characteristic function is K{t,n) = 
\og(t){t,n) = nlog(/o(2te)) -I- tz. The saddlepoint is the 
point t{v,n) that satisfies 



K'{t, n) 



2neli{2te) 



+ z 



(23) 



Eq. ([25)) is implicit and cannot be solved analytically. To 
derive an explicit equation we expand K'{t,n) in powers 
of e to second order. Doing so we get. 



t{v, n) 



2ne2 



(24) 



Substituting Eq. ([24]) into Eq. ([12]) we have. 



p{v,n\z,Q) 



1 ^(iz 

— e 2r 

2e 



-, 1/2 



(25) 



In Fig. [3| we present the exact numerical solution of given by Eq. ([^5]) [blue solid line] for n = {3,5,10,31}, 
Eq. ([5]) [red solid line] as well as the approximate one using only the second branch of the one-step TPF 
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FIG. 3. Histogram — upright crosses — of particle velocities 
after n — {3, 5, 10, 31} collisions, obtained by the iteration of 
Eqs. p3p . on the basis of an ensemble of 1.2 x lO" particles 
initially distributed as p{v,0) — S{v — e). The exact numerical 
solution of Eq. ((8]) [red solid line] as well as the approximate 
one given by Eq. pS)) [blue solid line] for n = {3, 5, 10, 31}, 
using only the second branch of the one-step TPF [Eq. (|19p ]. 
followed by the application of the method of images [Eq. (|20|l ] 
are also plotted for the sake of comparison, line. 



[Eq. (|T9l) ]. followed by the application of the method 
of images [Eq. (I^Ul) ]. The numerical solution is in total 
agreement with the histogram of particle velocities ob- 
tained by the iteration of the dynamical equations (jl3p 



[upright crosses], for all times. Even more, we see that the 
saddlepoint approximative solution describes very accu- 
rately the evolution of the PDF for n > 5. As can be ob- 
served, the PDF of particle velocities quickly approaches 
to a Gaussian distribution, in accordance with the pre- 
diction of the CLT. This is attributed to the fact that the 
TPF can be reduced to a difference kernel. Consequently, 
the additional assumption we made for the application of 
the CLT in Sec. lIII All namely that the statistical weight 
of the rare events^ is negUgible, is redundant. This can 
be circumvented, as aforementioned, by extending the 
domain of particle velocities. Thus, if one applies the 
CLT on the whole real line, then all the conditions for 
its applications are met exactly. As a final remark, we 
would like to stress that the success of the Fokker-Planck 
type of equation reported in the literature for even short 
times is attributed to the validity of the CLT, guarantee- 
ing that the PDF will converge to a normal distribution, 
allowing for the use of a diffusion equation. If however, 
the reduction of the TPF to a difference kernel is not fea- 
sible, then the transient can be arbitrarily long, a point 
demonstrated via an example in the following section. 



C. Long Transients 

In the last section we showed that the specific choice 
made for treating negative velocities after a collision, i.e. 
reflection with respect to the v = 0, reduces the TPF to 
an even function of the jump size. As a consequence, the 
PDF of particle velocities approaches rapidly to a sum 
of two spreading Gaussians. Clearly, after a number of 
collisions the system will "forget" its initial distribution, 
and the sum will converge to a single half-Gaussian cen- 
tered at w = 0+. Therefore, the most probable velocity 
for a particle evolving in the phase-randomized SFUM 
will eventually be Vp = 0"^, in clear contrast with the 
results given by the numerical simulation and analytical 
results derived using the exact dynamical mappin g^^''^^ , 
according to which as I) — >■ 0, p{v, n) — 0. From a phys- 
ical point of view this happens because if the motion of 
the wall in the configuration space is taken into account, 
as w — > collisions resulting in an energy loss can occur 
only in a small neighborhood around the wall's extreme 
positions, where its velocity is zero, resulting to a min- 
imal energy loss. Furthermore, if the particle velocity 
is comparable to the wall velocity, consecutive collisions 
can take place, resulting in a higher exit velocity from 
the interaction region within the exact model. 

On account of these properties of the collision process 
in the exact model, the refiection of negative velocities is 
not realistic. To gap the difference between the results of 
the simplified and the exact FUM, we propose instead of 
the inversion of negative particle velocities, the inversion 
of the direction of the wall's velocity, if the collision would 
lead to a negative particle velocity. This would lead in a 
greater energy gain in comparison with the refiection, as 
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ji; + u| < + Therefore, Eqs (|13p change to, 

2 



tn-1 + 

K-1 +2|u„| 
ecos(i„ + rin). 



(26a) 

(26b) 
(26c) 



Let us now derive the TPF. From Eqs. dTTJ), (US]) and 
(|26b[) we obtam, 



W{v, z) = -^pji^=^{e(2e + z-v) + Q{e- z/2) [Q{z - v - 2e) + Q{v - 2z)] } 
7r^4e2 - (u - z)2 ' 



(27) 



A comparison of the analytical result given by Eq. (j27p 
and the histogram of velocities obtained on the basis of 
an ensemble of 1.2 x 10^ particles after 1 iteration of 
Eqs. ([26]) is presented in Fig HI proving the validity of 
the derived result. 

As expected, the TPF has two branches, one taking 
effect only in the low velocity regime, i.e. z < 2e and 
another which is nonzero for any velocity z prior to a 
collision. As was also the case in the SFUM with re- 
flection of negative velocities, the part of the TPF that 
is relevant to the low-velocity regime, depends on the 
jump size, as well as, on the velocity of a particle prior 
to a collision. However, due to the fact that this branch 
of W does not have a simple geometrical interpretation, 
the single-step transition function cannot be reduced to 
a difference propagator by an extension of the domain of 
W to the whole real line. Thus, the conditions for the 
application of the CLT are not met exactly. However, 
due to the acceleration of the particles, as n — oo, the 
probability measure of particles having velocity z < 2e 
becomes negligible. Therefore, for n ^ 1 and t; ^ e the 
PDF of particle velocities tends to a Gaussian distribu- 
tion [Eq. (dH)]. 

The study of the transient behaviour of the PDF re- 
quires the solution of the CKE [Eq. The numerical 
solution of Eq. ([S]) at times n = {3,5,17,316} is pre- 
sented in Fig. [5] The histograms of particle velocities for 
the same times, calculated by iterating an ensemble of 
1.2 X 10^ particles for up to n = 10^ collisions, are also 
plotted for the sake of comparison. It can be seen, that 
the solution of the CKE is in agreement with the results 
of the simulation for all times presented. In Fig. [HI the 
histogram of velocities for n = 10^ collisions is plotted. 
The solution obtained from the application of the CLT 
— on the assumption that the statistical weight of col- 
lisions happening in the region v < 2e is negligible, is 
also plotted, and is full agreement with the PDF in this 
velocity region. However, a blow-up of the low-velocity 
region shows that even after 10^ collisions, the PDF di- 
verges from the Gaussian profile. This is clear evidence 
that even after a very large number of collisions, the PDF 
in the whole velocity domain cannot be described by an 




0.1 0.2 0.3 0.4 



FIG. 4. Histogram — diagonal crosses — of particle velocities 
after a single collision with the "moving" wall, obtained using 
Eqs. (|26|) and an ensemble of 1.2 x 10^ particles with initial 
velocity z — 0.1. The analytical result [Eq. (|27|l ] for the 
one-step transition probability is also plotted for the sake of 
comparison — solid line 



EPF, in contrast to the standard version of the SFUM 
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FIG. 6. Histogram — upright crosses — of particle velocities 
after n = 10^ collisions, obtained by the iteration of Eqs. ((26}, 
on the basis of an ensemble of 1.2 x 10® particles initially dis- 
tributed as p(v,0) = S{v — e). The asymptotic Gaussian dis- 
tribution [Eq. (|14|) ] predicted by the CLT is also plotted -solid 
(red) line. A blow-up of the numerically obtained histogram 
at the low- velocity region is illustrated in the inset. 



D. Fermi acceleration in the exact FUlVI 



FIG. 5. Histogram — upright crosses — of particle velocities 
after n = {3, 5, 17, 316} collisions, obtained by the iteration 
of Eqs. (|26p . on the basis of an ensemble of 1.2 x 10® particles 
initially distributed as p{v,0) = S{v — e). The solution ob- 
tained by numerically solving the forward CKE [Eq[8] using 
as the TPF of the modified version of the SFUM [Eq. [27] is 
also plotted — solid (red) line. 



[Eqs. ([T^] even though the argumentation used in both 
cases was the same, i.e. the particles are accelerated. 
This exemplifies the potential pitfalls of a diffusion ap- 
proximation of Fermi-acceleration in time-dependent bil- 
Uards. 



The simplification employed to the treatment of Fermi 
acceleration in the FUM — treating the wall as fixed in 
real space — however widespread and thoroughly studied, 
prohibits the study of the details of Fermi acceleration. 
In RefSfi^i^'^ it was proved that small additional fluctu- 
ations of the time of collision due to dynamical correla- 
tions induced by the displacement of the scatterer upon 
impact, quantitatively as well as qualitatively change the 
evolution of the PDF of velocities, increasing the effi- 
ciency of Fermi acceleration. Moreover, the development 
of correlations causes the CLT to break down and the 
asymptotic PDF ceases to be a normal distribution. 

In Ref.^'* utilizing a novel simplification, the so-called 
hopping approximation, which succeeds into retaining all 
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the essentials of the exact dynamics, an analytical solu- 
tion for the asymptotic behaviour of the PDF of particle 
velocities in the exact FUM, which was in excellent agree- 
ment with the numerical simulation of the exact FUM. 
Specifically, it was shown by means of a Fokker-Planck 
type of equation that, in contrast with the SFUM, the 
attractor of the PDF of velocities in the function space is 
a Maxwell-Boltzman like distribution, i.e. independently 
of the initial distribution of velocities, the PDF converges 
to a Maxwell-Boltzman like distribution. Therefore, in 
the case of the exact FUM for v — ?> 0+, p{v, n\z, 0) — > 0, 
in contrast to the SFUM where p{v,n\z,0) attains its 
maximum value for w — >■ 0"*". This difference between the 
simplified and the exact FUM can be understood as fol- 
lows: If the velocity of a particle after a collision with the 
moving wall is small, then multiple successive collisions 
are likely to occur within the exact FUM, resulting into 
higher exit velocities, as opposed to the simplified model, 
within which successive collisions cannot be realized. 

Another subtle difference between the simplified and 
the exact FUM [see Sec. IIII A 3| is that within the exact 
FUM particles with low velocity are more likely to collide 
with the oscillating wall near its turning points, where 
the wall velocity is close to 0. Hence, the velocity jump 
performed by a particle due to a collision with the wall 
Av — as w — 0. As a result, Fermi acceleration when 
using the exact dynamics can be better approximated by 
a continuous stochastic process, or equivalently, by the 
FPE. Still, the transient statistics in the system can only 
be studied by means of the CKE. 

As aforementioned, the movement of the wall in the 
configuration space described by the exact dynamics re- 
sults into a more efficient energy transfer from the mov- 
ing wall to the particles upon collision, when compared to 
the SFUM. In mathematical terms, this causes the PDF 
of the oscillation phase on collisions to deviate from the 
uniform distribution, reffecting the fact that head-on col- 
lisions are more preferable than head-tail collisions. How- 
ever, the phase of oscillation of the moving wall when a 
particles collides with the fixed wall — or when it passes 
through any fixed point within the area between the two 
walls comprising the FUM — is uniformly distributed. 
The map describing the exact dynamics is, 

dn, = e sin {St„ + i„-i + (28a) 
Un = e cos{Stn + <„_i + rjn) (28b) 

Vn = Vn-1 + 2Un, (28c) 

where d„. stands for the position of the moving wall in 
the instant of the nth collision, w„ for the wall velocity, ry„ 
for the random phase component and w„ for the particle 
velocity after the 71th collision. The time of free flight 
6tn is obtained by solving the implicit equation 

Xn-l + Vn-lSt„ = dn, (29) 

where a;„ stands for the position of the particle in the 
instant of the nth collision. 
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FIG. 7. Histogram — upright crosses — of particle velocities 
after n = {3, 5, 17, 316} collisions, obtained by the iteration 
of Eqs. (|28p . on the basis of an ensemble of 1.2 x 10® particles 
initially distributed as p{v,0) — 5{v — 10). The solution ob- 
tained by numerically solving the forward CKE [Eq|8] using 
as the TPF of the exact model [Eq. [31] is also plotted — solid 
(red) line. 

If we denote the phase of oscillation of the moving wall 
when a particle collides with the fixed wall with -0, then 

= cos-1 (^) + ^ (l + V^^^) (30) 
For Eq. (1501) we obtain for the distribution of the wall 
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velocity upon collision, 



Pe{u) 



u + z 



(31) 



Installing Eq. (PT|) into Eq. (ITS)) we obtain the one-step 
TPF for the exact model 



Q(2t-\v-z\){v^^z) 



27rzy4e 



zf 



(32) 



Inserting Eq. ([5^ into the CKE we can numerically com- 
pute the evolution of the PDF of particle velocities. In 
Fig. [3 the numerical solution of the CKE is compared 
with the histogram of particle velocities, obtained by 
simulating 1.2 x 10^ trajectories using Eqs. (|28p . The 
particles were initially distributed according to p{v, 0) = 
8{y— 10). Once more, the solution of the CKE is in com- 
plete agreement with the results of the simulation, prov- 
ing that the method can also be successfully employed 
when the exact dynamics are taken into account. 



IV. SUMMARY AND CONCLUSIONS 

Fermi acceleration, is one of the most interesting as- 
pects of time-dependent billiards, as it has been under- 
stood over the years, that it is a fundamental acceleration 
mechanism, playing a key role in a variety of phenomena, 
far beyond its original scope, i.e. cosmic-ray particle ac- 
celeration. 

Until now, the investigation of Fermi acceleration 
— in the class of time-dependent billiards in which it 
develops — has been carried out via its approximation 
with a diffusion process. Within, this framework, the evo- 
lution of the density of particle velocities was determined 
by a Fokker-Planck equation. However, its derivation is 
always based on assumptions and approximations that 
rarely can be justified. Moreover, its prediction power is 
limited in the long-term statistics of the system and no 
information is given for the transient behaviour. Even 
more, its use in the SFUM, which is the first system that 
was successfully investigated with the use of the FPE, is 
completely redundant, as the CLT yields the same results 
in a far more straightforward manner. 

Herein, we proposed a consistent methodology, which 
obviates unclear assumptions and even more, can give 
an accurate description of the transient evolution of par- 
ticle velocities. The cornerstone of this methodology is 
the use of the Chapman-Kolmogorov equation. The fun- 
damental difference in comparison with the traditional 
approach using the FPE, is that no assumption for the 
continuity of the stochastic process describing Fermi ac- 
celeration needs not to be made. Another advantage of 
the proposed approach is that all collision events can be 
taken into account, which cannot be done in the construc- 
tion of the FPE, and even when possible, it can lead to 
less accurate results. 



The method was successfully applied to the FUM, 
which is the prototype of time-dependent billiards ex- 
hibiting Fermi acceleration. In specific, we studied the 
standard SFUM, within which collisions leading to a po- 
tential escape from the system are handled by artificially 
inverting the particle exit velocity, as well as a variant 
of the SFUM, where if a collision would lead to a par- 
ticle still moving towards the wall, the velocity of the 
wall is inverted before the collision takes place. Finally, 
we showed how this method can be applied to the exact 
model, showing how the effect of the motion of the wall 
in the configuration space can be included in the descrip- 
tion of Fermi acceleration through the CKE. In all three 
cases, the CKE yielded accurate results for all times. As 
a final remark, we would like to stress that this methodol- 
ogy can be applied to higher-dimensional billiards^ and 
therefore is generic. 
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